      SUBROUTINE SETUP
      RETURN
      END
C   HS73
      BLOCK DATA
      INCLUDE 'BLOCO.INC'
      INTEGER I,J
      DATA NAME/'HS73'/
      DATA (X(I),I=1,4)/4*1.D0/
      DATA N/ 4/ , NH/1/ , NG/ 6/
      DATA DEL0/1.00D0/ ,TAU0/1.D-1/ ,TAU/.1D0/
      DATA (GUNIT(1,J),J=4,7)/4*1/,(GUNIT(3,I),I=4,7)/4*1/,
     F     (GUNIT(2,I),I=4,7)/1,2,3,4/
     F     ,((GUNIT(I,J),I=1,3),J=0,3)/-1,0,0,-1,0,0,-1,0,0,-1,0,0/
      END
      SUBROUTINE EF(X,FX)
      INCLUDE 'FUCO.INC'
      INTEGER J
      DOUBLE PRECISION X(NX),A(4),FX,S
      SAVE A
      DATA  A/24.55D0,26.75D0,39.0D0,40.5D0/
      ICF=ICF+1
      S=0.D0
      DO      100      J=1,4
      S=S+A(J)*X(J)
  100 CONTINUE
      FX=S
      RETURN
      END
      SUBROUTINE EGRADF(X,GRADF)
      INCLUDE 'FUCO.INC'
      INTEGER J
      DOUBLE PRECISION X(NX),GRADF(NX),A(4)
      SAVE A
      DATA  A/24.55D0,26.75D0,39.0D0,40.5D0/
      ICGF=ICGF+1
      DO      100      J=1,4
      GRADF(J)= A(J)
 100  CONTINUE
      RETURN
      END
      SUBROUTINE EH(I,X,HXI)
      INCLUDE 'FUCO.INC'
      INTEGER I
      DOUBLE PRECISION X(NX),HXI
      CRES(I)=CRES(I)+1
      HXI=X(1)+X(2)+X(3)+X(4)-1.D0
      RETURN
      END
      SUBROUTINE EGRADH(I,X,GRADHI)
      INCLUDE 'FUCO.INC'
      INTEGER I,J
      DOUBLE PRECISION X(NX),GRADHI(NX)
      CGRES(I)=CGRES(I)+1
      DO      100       J=1,4
      GRADHI(J)=1.D0
  100 CONTINUE
      RETURN
      END
      SUBROUTINE EG(I,X,GXI)
      INCLUDE 'FUCO.INC'
      INTEGER I,J
      DOUBLE PRECISION X(NX),A(4),B(4),C(4),GXI,S,SL
      SAVE A,B,C
      DATA B/12.D0,11.9D0,41.8D0,52.1D0/
      DATA A/2.3D0,5.6D0,11.1D0,1.3D0/
      DATA C/.28D0,.19D0,20.5D0,.62D0/
      IF ( GUNIT(1,I+NH) .EQ. -1 ) CRES(I+NH)=CRES(I+NH)+1
      IF(I .GT. 2)      GOTO 400
      IF(I .EQ. 2)      GOTO 200
      S=-5.D0
      DO      100      J=1,4
      S=S+A(J)*X(J)
 100  CONTINUE
      GXI=S
      RETURN
  200 CONTINUE
      S=0.D0
      SL=-21.D0
      DO      300      J=1,4
      SL=SL+B(J)*X(J)
      S=S+C(J)*X(J)**2
  300 CONTINUE
      GXI=SL-1.645D0*DSQRT(S)
      RETURN
  400 CONTINUE
      GXI=X(I- 2)
      END
      SUBROUTINE EGRADG(I,X,GRADGI)
      INCLUDE 'FUCO.INC'
      INTEGER I,J
      DOUBLE PRECISION X(NX),GRADGI(NX),B(4),A(4),C(4),S
      SAVE A,B,C
      DATA B/12.D0,11.9D0,41.8D0,52.1D0/
      DATA A/2.3D0,5.6D0,11.1D0,1.3D0/
      DATA C/.28D0,.19D0,20.5D0,.62D0/
      IF(I .GT. 2)    GOTO 400
      IF(I .EQ. 2)    GOTO 200
      DO     100      J=1,4
      GRADGI(J)=A(J)
  100 CONTINUE
      RETURN
  200 CONTINUE
      S=0.D0
      DO      250      J=1,4
      S=S+C(J)*X(J)**2
  250 CONTINUE
      S=-1.645D0/DSQRT(S)
      DO      300      J=1,4
      GRADGI(J)=B(J)+S*C(J)*X(J)
  300 CONTINUE
      RETURN
  400 CONTINUE
      DO      500    J=1,4
      GRADGI(J)=0.D0
  500 CONTINUE
      GRADGI(I-2)=1.D0
      RETURN
      END


